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Abstract 


A  signal  processing  technique  is  presented  for  correcting  imbalances  and 
distortions  introduced  to  the  signals  by  the  phase  detectors  of  a  coher¬ 
ent  wideband  radar.  The  signal  model  and  sources  of  signal  errors  are  de¬ 
scribed,  an  analytic  description  of  how  the  corrections  are  derived  is  pro¬ 
vided,  and  a  sample  application  is  presented  with  the  use  of  simulated 
data.  This  report  has  been  prepared  in  anticipation  of  a  subsequent  re¬ 
port  in  which  the  performance  of  this  signal  processing  technique  will  be 
compared  with  the  performances  of  several  other  techniques  developed  for 
similar  purposes  with  actual  data  measured  in  the  field. 
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1.  Introduction 


The  U.S.  Army  Research  Laboratory  (ARL)  has  a  variety  of  instrumenta¬ 
tion  radars  for  research  related  to  sensor  technologies,  to  include  model¬ 
ing,  simulation,  analysis,  and  signal  processing  of  radar  signatures.  Before 
any  of  these  research  objectives  are  met,  the  data  collected  by  the  instru¬ 
mentation  radars  must  be  calibrated  to  relate  the  returned  signal  to  the 
transmitted  signal.  This  calibration  removes  errors  introduced  by  compo¬ 
nents  of  the  radar  and  scales  the  data.  In  particular,  ARL's  Millimeter- Wave 
Branch  of  the  Radio  Frequency  (rf)  and  Electronics  Division  has  several  in¬ 
verse  synthetic  aperture  radars  containing  phase  detectors  that  require  er¬ 
ror  correction.  As  these  new  radars  are  developed  and  introduced,  new  cal¬ 
ibration  techniques  have  been  developed  so  that  three  distinct  techniques 
are  now  being  used.  Two  of  these  techniques,  the  method  of  Wallace  and 
Pizzillo1  and  a  phase  modulation  and  demodulation  technique,2  have  been 
documented  as  ARL  technical  reports.  The  third,  a  harmonic  decomposi¬ 
tion  technique,  has  only  been  documented  as  an  algorithm  in  software. 

This  report  provides  an  analytic  description  of  the  third  algorithm  in  antici¬ 
pation  of  publishing  a  report  detailing  a  comparison  of  the  three  techniques 
with  the  use  of  data  collected  by  the  state-of-the-art  35-GFIz  monopulse  in¬ 
strumentation  radar.  This  instrument  uses  the  phase  modulation  and  de¬ 
modulation  technique,  an  architecture-based  calibration,  and  the  only  one 
of  the  various  radars  that  can  collect  data  that  may  be  corrected  by  all  three 
techniques,  thereby  allowing  a  direct  comparison  of  the  efficiency  and  ac¬ 
curacy  of  each  scheme. 


1H.  Bruce  Wallace  and  Thomas  J.  Pizzillo,  A  Technique  for  Calibrating  the  Phase  Detector 
of  a  Wideband  Radar  Using  an  External  Target,  U.S.  Army  Research  Laboratory,  ARL-TR-1521 
(March  1998). 

2Thomas  J.  Pizzillo  and  H.  Bruce  Wallace,  A  Technique  for  Calibrating  the  Phase  Detector 
of  Wideband  Radars  Using  a  Phase  Modulation  and  Demodulation  Scheme,  U.S.  Army  Research 
Laboratory,  ARL-TR-1567  (May  1998). 
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2.  Derivation  of  Correction  Terms 


Figure  1  is  a  simplified  block  diagram  of  a  polarimetric,  two-coordinate, 
amplitude-comparison  monopulse  system  operating  at  94  GHz,  with  a 
640-MHz  bandwidth.  Either  linearly  or  circularly  polarized  radiation  is 
transmitted  and  both  components  are  received,  i.e.,  left-circular  transmit 
and  right-  and  left-circular  receive  or  vertical  transmit  and  vertical  and 
horizontal  receive.  Transmit  polarization  may  be  changed  to  pulse  to  pulse 
or  ramp  to  ramp. 

The  rf  local  oscillator  (LO)  is  frequency-stepped  synchronously  with  the 
transmitter  and  maintains  a  constant  3-GHz  offset.  This  source  is  mixed 
with  the  received  signal  to  provide  the  3-GHz  intermediate  frequency  (IF) 
to  the  in-phase/ quadrature  (1/ Q)  detectors.  The  IF  is  mixed  with  the  3-GHz 
LO  to  provide  the  final  direct  current  (DC)  signal  to  the  analog-to-digital 
converter  (ADC).  The  3-GHz  LO  is  also  the  source  of  the  injected  test  signal 
used  for  correcting  the  phase  and  gain  imbalances  introduced  by  the  I/Q 
detectors.  A  simplified  block  diagram  of  the  IF  section  is  shown  in  figure  2. 
This  test  signal  is 

S  =  Acos27t/  ,  (1) 

where  /  represents  the  test  signal's  constant  frequency  and  A  is  the  am¬ 
plitude  of  the  test  signal.  However,  this  only  provides  a  single  DC  output 


Figure  1.  Simplified 
block  diagram  of  a 
monopulse  radar. 
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3-GHz  test  signal 


Figure  2.  Simplified 
block  diagram  of  IF 
section  of  a  monopulse 
radar. 


value  for  a  given  frequency  from  the  I/Q  detectors,  and  correcting  the  out¬ 
put  of  the  detectors  over  a  range  of  frequencies  is  necessary.  To  accomplish 
this,  one  must  phase-modulate  the  test  signal  to  provide  an  input  to  the 
I/Q  detectors  as 

S(m)  =  A  cos  (2-Kf +  2TT^j  ,  (2) 

where  M  is  an  even  integer  representing  the  number  of  phase  modula¬ 
tion  steps.  Shifting  the  phase  creates  a  pseudo-time-sampled  signal  with 
frequency  of  1  Hz  and  a  period  of  1  /  M.  This  signal  is  injected  at  the  IF 
section  and  then  propagates  through  the  phase  detector,  where  a  second 
signal,  shifted  90°  relative  to  the  test  signal,  is  generated.  These  signals  are 
referred  to  as  the  I  and  Q  channels  of  the  detector: 

I(m)  =  A  cos  +  Vi  and 

Q(m)  =  CM  sin  +  <5^  +  Vq  ,  (3) 

where  V  and  Vq  are  DC-offsets  and  G  and  8  represent  the  relative  gain  and 
phase  imbalances  between  the  outputs  of  I  and  Q  channels,  all  introduced 
by  the  imperfect  detector.  The  2irf  term  has  been  eliminated  because  the 
defection  is  a  mixing  process  with  /  as  the  baseband.  At  this  juncture,  these 
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errors  could  be  corrected  by  the  method  described  in  Wallace  and  Pizzillo. 
However,  in  addition  to  the  errors  introduced  by  the  phase  detector,  the 
components  that  compose  the  IF  section  cause  nonlinear  distortion  of  the 
signal  so  that  the  signals  may  be  better  modeled  by  a  set  of  trigonometric 
polynomials, 

I(M)  =  ]T  A  cos  +Vt  N  =  0.  2. 4 . . .  ^ 

r?=o  \  J\1  /  2 

Q(^)  =  B”  sin  (2n7TJj)  +  VQ  N  =  1-3.5,...  (¥■  -  l)  .  (4) 

Although  the  phase  and  gain  imbalance  of  the  phase  detector  creates  a  re¬ 
sponse  at  the  image  range  bin,3  this  error  is  subsumed  by  the  appropriate 
terms  and  coefficients  of  the  polynomial  and  needs  not  be  explicitly  enu¬ 
merated.  Equation  (4)  represents  a  set  of  functions  orthogonal  on  the  com¬ 
plex  plane  over  the  range  [— 7r,  n].  Hence,  the  coefficients  A  and  Bn  may  be 
determined  via  the  complex  fast  Fourier  transform  (CFFT)  of  the  measured 
1/ Q  detector  outputs  by 

F  [I  (M) ,  Q  (A/)]  =  |  A  cos  ^2nn — ^  +  iBn  sin  ^n7r— ^  +  V  j  exp 

=  I  (K)  +  iQ  (K) ,  (5) 

where  V  =  V,  +  Vq  and  I (K)  and  Q (K)  are  the  M  outputs  of  the  CFFT  that 
represent  the  harmonic  coefficients  of  the  transformed  signal.  That  is,  for 
each  k,  only  the  trigonometric  component  k  —  n  contributes 

I  (K)  =  A  +  A  cos  ^27t— ^  +  . . .  +  Am -2  cos  ,  and 

Q  (K)  =  B\  sin  +  S3  sin  +  •  •  ■  +  BM^l  sin  ^ 

The  maximum  number  of  contributing  harmonics  is  determined  by  the 
Nyquist  rate  or  M/2  distributed  symmetrically  about  the  fundamental: 

I  (K)  =  Ao  H~  A cos  (— 7r)  +  . . .  +  Am.  cos  (tt)  ,  and 

2  2 

Q  (K)  =  B 2-_ki_  sin  ( 27r= 2^-)  +  •  •  •  +  Bhl^l  sin  •  (7) 

By  repeating  this  procedure  for  several  amplitudes  and  storing  the  values 
in  a  lookup  table,  one  may  correct  any  measured  amplitude  on  the  complex 
plane  defined  by  I  and  Q  up  to  the  maximum  calibration  amplitude. 

’Merril  Skolnik,  Radar  Handbook,  2nd  ed.,  McGraw-Hill,  Inc.  (1990),  p  3.41. 
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3.  Technique  Application 


Figure  3.  Simplified 
block  diagram  of  rf 
section  of  a  monopulse 
radar. 


The  correction  noted  in  section  2  is  applied  in  the  following  procedure:  A 
constant  amplitude  signal  is  injected  into  the  IF  section  of  the  radar  via  the 
toggle  switches  (see  fig.  2).  This  provides  a  DC  signal  to  the  ADC.  Data 
are  collected  at  16  phase  angles,  22.5°m,  where  m  =  0,1. . .  15,  controlled  by 
the  phase-shifter  shown  in  figure  3  for  both  the  I  and  Q  channels.  Six  thou¬ 
sand  four  hundred  samples  are  averaged  for  each  phase  setting  to  mini¬ 
mize  noise.  The  data  are  converted  to  signed  integer  values  by  subtracting 
half  the  dynamic  range  of  the  ADC  from  each  set,  i.e.,  2048.  A  complex 
value  is  generated  from  the  I  and  Q  data  and  then  sorted  from  the  mini¬ 
mum  to  maximum  phase.  This  process  is  repeated  for  four  attenuation  lev¬ 
els:  98,  78,  58,  and  39  percent  of  the  A/D  dynamic  range,  i.e.,  4095.  The 
final  results  are  four  sets  of  complex  data  of  16  values  each  from  which  the 
correction  coefficients  are  derived. 

The  correction  coefficients  are  used  to  correct  measured  data  in  both  phase 
and  amplitude  by  subtracting  the  unwanted  harmonic  components  from 
each  measured  data  point,  Jm  and  Qm.  Measured  amplitudes  are  corrected 
with  the  correction  amplitudes  above  and  below  the  measured  value,  with 


94-GHz  transmit  signal 
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an  appropriate  weighting  factor.  This  weighting  factor  is  calculated  with 


ur  _  (Ay  Am)  TJr  (Am  —  A[) 

““  (A,  -A,)  '  “  (A,  -  A)  '■ 


(Wv  +  w,  =  1)  . 


(8) 


where  Au  is  the  correction  amplitude  greater  than  the  measured  value. 
At  is  the  correction  amplitude  less  than  the  measured  value,  and  Am  - 
A;n  =  the  amplitude  of  the  measured  I  and  Q  data  points. 

Consider  the  simulated  data  in  figure  4  created  for  M  =  16.  The  set  of  meas¬ 
ured  I/Q  pairs  is  indicated  with  squares  and  falls  between  the  calibration 
amplitudes,  Au  =  1500  and  Aj  =  1200.  The  corrected  I  and  Q  values  are 
determined  by 


Ic  —  Im  W' u  [-dou  +  A_8tt  COS  +  .  .  .  +  Aft,,  cos  (86 m )] 

-  W\  [Aqi  +  t4_8/  cos  (—8 0m)  +  . . .  +  Agi  cos  (80,„)]  .  and 
Qc  —  Qm  ~  Wv  [^4-7u  sin  (—7dm)  +  . . .  +  Ajv  sin  (7 Qm)\ 

-  Wi  [A_7/  sin  (-7  0m)  +  ...  +  Aji  sin  (7Qm)] .  (9) 

where  Ic  and  Qc  are  the  corrected  I  and  Q  data  points.  6m  =  arctan  2  (j1^ 
is  the  measured  relative  phase  between  Im  and  Qrn.  Anu  and  An,  are  the  nth 
harmonic  coefficient  for  the  upper  and  lower  calibration  bounds,  respec¬ 
tively.  These  expressions  are  applied  to  each  of  the  16  I/Q  pairs  indicated 
with  circles  in  figure  4.  This  example  is  for  demonstration  only  and  for  clari¬ 
fying  the  application  of  this  technique.  A  subsequent  report  will  compare 
the  application  of  this  technique  to  measured  data,  and  its  performance 
will  be  compared  with  two  other  techniques  developed  for  the  same  data 
correction  problem. 


Figure  4.  Simulated  I/Q 
data:  Squares  indicate 
measured  data  and 
circles  indicate 
corrected  data. 
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4.  Conclusion 


An  algorithm  currently  used  by  ARL  (Millimeter- Wave  Branch)  for  correct¬ 
ing  the  I  and  Q  errors  introduced  by  the  phase  detectors  of  a  wideband 
radar  has  been  presented.  The  signals  with  errors  have  been  modeled  as 
trigonometric  polynomials  and  the  correction  described  as  a  harmonic  de¬ 
composition  based  upon  the  CFFT.  An  analysis  of  the  performance  of  this 
technique  is  left  for  a  subsequent  publication  in  which  the  performance  of 
the  two  other  techniques  (see  sect.  1)  will  be  compared  along  with  the  re¬ 
sults  of  this  harmonic  decomposition  as  applied  to  actual  measured  data 
from  a  state-of-the-art  35-GHz  monopulse  instrumentation  radar. 
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